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^ ' Abstract 

■ We extend the theory of leading twist nuclear shadowing to calculate leading twist nuclear 
K*" . difFractive parton distribution functions (nDPDFs). We observe that the quark and gluon nPDFs 

\ have different patterns of the ^-dependence. It is found that the probability of diffraction in the 

quark channel increases with A, reaching about 30% at x ~ 10~^ for A ~ 200, and weakly decreases 

■ with Q^. In the gluon channel, the probability of diffraction is large for all nuclei (~ 40% for heavy 
' nuclei at x 10""^ and Qq ~ 4 GeV^), it weakly depends on A and it decreases rather fast with 
. increasing - the probability decreases by approximately a factor of two as changes from 
\ 4 GeV'^ to 100 GeV'^. We also find that nuclear shadowing breaks down Regge factorization of 

nDPDFs, which is satisfied experimentally in the nucleon case. All these novel effects in nDPDFs 



Oh! 



are large enough to be straightforwardly measured in ultraperipheral collisions at the LHC. 



1 Introduction 

I It is firmly established by the HERA HI |2] and ZEUS [3] experiments that inclusive 
diffraction constitutes a significant fraction, about 10%, of the total cross section for deep 
inelastic scattering (DIS) of leptons on hydrogen. The diffractive events are characterized by 
the absence of hadronic activity in the detector at central rapidities. This QCD phenomenon 
involves both nonperturbative and perturbative aspects of the QCD dynamics, for the dis- 
cussions see and it is usually referred to as scattering off the "Pomeron". Note that in 
the case of DIS, the dynamics is quite different from that of the Pomeron pole exchange in 
soft hadron-hadron interactions, so we use the concept of Pomeron in order to indicate that 
we refer to the kinematics of large rapidity gaps and small momentum transfers to the recoil 
nucleon. 

The proof of the factorization theorem for hard diffraction |2I enables one to describe 
the process in terms of Q^-dependent diffractive parton distribution functions (DPDFs) and 
extract the DPDFs from various diffractive data. The current data 13 El are consistent 
with the dominance of leading twist in hard diffraction and with the dominance of the gluon 
DPDF over the sum of the quark DPDFs. The DPDFs for the scattering off unpolarized 
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target, fj/j^ {(3, Q^, xjp, t) depend on four variables: Bjorken x, virtuality Q^, four-momentum 
transfer squared to the target t and the fraction of the longitudinal momentum loss by the 
target xjp. 

The aim of the present paper is to investigate, within the leading twist approximation, nu- 
clear dependence of coherent (without nuclear break-up) diffraction induced by hard probes 
and to obtain nuclear DPDFs. This would allow to calculate the cross sections of various 
diffractive DIS processes as well as of direct photon diffraction off nuclear targets. Experi- 
mental studies of these processes will be feasible in ultraperipheral collisions (UPC) of heavy 
ions at the LHC P and at the EIC 0. 

Inclusive hard diffraction in DIS off nuclei has been studied in a number of papers, see 
e.g. [HIE]- Unlike all previous attempts, we use the QCD factorization theorem for hard 
diffraction jH] and the leading twist theory of nuclear shadowing fUl ^] • This enables 
us, for the first time, to calculate nDPDFs, i.e. to perform the flavor separation. This is an 
essential ingredient for the calculation of various, more complicated, problems such as charm 
production in DIS and dijet production in direct photon diffraction off nuclear targets, which 
need nDPDFs and especially the gluon nDPDF. The present paper complements our studies 
of the role of leading twist nuclear shadowing at small-x in inclusive processes with nuclei. 

While we assume that the QCD factorization theorem for hard diffraction in DIS holds 
and all considered effects are leading twist effects, the studies of coherent diffraction on nuclei 
help to understand the transition to the regime of high parton densities. Indeed, the fraction 
of diffraction of the total cross section in DIS is a measure of how close to the black body 
regime (the regime of complete absorption of the projectile by the target) one is. While 
diffraction is approximately 10% of the total cross section in DIS on hydrogen, the fraction 
of diffractive events steadily increases as one increases the atomic number A, asymptotically 
approaching the absolute limit of one half. 

This paper is structured as follows. In Sect. |21 we recapitulate essential points of the 
leading twist nuclear shadowing model and present the formula for coherent diffraction on 
nuclei. The analysis and discussion of the resulting expressions for nuclear diffractive parton 
distributions are presented in Sect. El We conclude and summarize in Sect. El 

2 Leading twist nuclear shadowing and coherent diffrac- 
tion on nuclei 

The theory of leading twist nuclear shadowing is based on the Gribov's relation between 
nuclear shadowing and diffraction |14 , Collins factorization theorem for hard diffraction in 
DIS j3] and the QCD analysis of the HERA data on hard diffraction in DIS on hydrogen 
|2|. The foundations of the resulting theory and predictions for nuclear parton distribution 
functions and inclusive structure functions can be found in Refs. [TIH IIH IT^. 

The master equation for the evaluation of the shadowing correction, Sfj/A, to the nuclear 
structure parton distribution functions of flavor j, fj/A = ^fj/N — ^fj/A has the form 
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with 7] the ratio of the real to imaginary parts of the diffractive scattering amphtude; zi, 
Z2 and h the longitudinal (in the direction of the incoming virtual photon) and transverse 
coordinates of the nucleons involved (defined with respect to the nuclear center); P, xjp 
and t the usual kinematic variables used in diffraction; /3 = x/xjp; tmin ~ 0; pA{b,Zi) the 
nucleon distribution in the target nucleus. The upper limit of integration, xjp^q is a cut-off 
parameter, which equals 0.1 for quarks and 0.03 for gluons. The effective cross section, 0"^^, 
is expressed through the nucleon DPDFs as (see Ref. [12] for the detailed discussion and 
numerical estimates) 

Equation (jT)) serves to define the input nuclear PDFs at the initial scale Ql, Ql = 4 GeV^ in 
our analysis. Nuclear PDFs at larger scales are obtained using the NLO QCD evolution 
equations. 

In Eq. (^, the interaction with two nucleons is calculated in a model-independent way. 
The only source of mo del- dependence is due to the approximation of the interaction with 
three and more nucleons by the attenuation factor e '■'^^^-"■"'^ *'')°"off dzpA(b,z) ^ which involves 
cXgfj, the rescattering cross section given by Eq. (0). While this quasi-eikonal approximation 
is expected to be valid ai Ql = A GeV^, it becomes progressively worse with increasing 
Q^. The reason for this is that the eikonal approximation conserves the number of bare 
particles and thus contradicts QCD evolution. As a result, one obtains a wrong, higher 
twist, Q^-dependence of nuclear shadowing in the processes dominated by small partonic 
configurations of the incoming virtual photon. Only at low scales, where the effects of 
QCD evolution are not very important, can one justify the use of the eikonal and quasi- 
eikonal approximations. This means that Eq. (Q) should be used only at the initial scale 
Ql = A GeV2. 

The generalization to the case of coherent diffraction in DIS on nuclei is rather straight- 
forward, and it follows closely the case of the vector meson diffraction, see e.g. The 
nuclear diffractive parton distribution of flavor j can be presented in the form 

f^/fi^^Ql^^ip) = ^167r/jy}J^(a;,Qo>a;jp,tmin) j d% 

2 

The superscripts (3) and (4) denote the dependence on three and four variables, respectively. 
We present our Eq. Q for the t-integrated nuclear DPDFs since it is more compact and since 
it is not feasible to measure t in diffraction off nuclei in the collider experiments. In deriving 
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Eq. we neglected a possible /3-dependence of cTgg(x, Q^) in the exponential factor and 
substituted a^^ by its average value. Since the total probability of diffraction changes rather 
weakly with Cgg, see e.g. jTHj, this seems a reasonable first approximation. At the same 
time, in the region of small (3 and small x corresponding to the triple Pomeron kinematics 
for soft inelastic diffraction, we expect a significant suppression of diffraction as compared 
to the quasi-eikonal approximation of Eq. Q for ^ Ql, see the discussion in the end of 
the section. 

One should note that the large momentum transfer Q^, which is necessary for the ap- 
plicability of the QCD factorization theorem, does not preclude the existence of coherent 
nuclear diffraction. Indeed, at high energies, the minimal momentum tranfer to the nucleus 
tmin is small, tinin ~ x'^jM^, which makes it possible for nucleus to stay intact (or diffract 
into low mass excited states). In practice, coherent nuclear diffraction can be identified by 
its distinctly sharp t-dependence in the forward direction (forward diffractive peak), which 
originates from the factor (F4(t))^ where FA_{t) is the nuclear form factor. 

Assuming the exponential t-dependence of inclusive diffraction on free nucleons and using 
that tmin ~ 0, we obtain 

/^JJ^x, Ql xjp, t = 0)= ff/^\x, Ql xjp) , (4) 

where Bj is a slope. In our analysis, we use Bj = 7.2 GeV~^ for quarks and Bj = 6 + 
0.25 ln(10~^/x) GeV~^ for gluons, which are rather close numerically. 

For sufficiently small values of xjp, xjp < 0.01, the HI data P3I2] and ZEUS data ^ can 
be fitted reasonably well using a factorized approximation first suggested within the picture 
of soft mechanism of diffractive process by Ingelman and Schlein jT7]. In this approach, 
DPDFs can be presented as a product of a factor depending only on t and xjp (Pomeron 
flux) and a factor depending only on [3 = x/xjp and Q^ (which is often referred to as the 
DPDF of the Pomeron) 

f^/N\^'Ql^^ip) = fp/pixip) fj/ip{l3 = x/xip,Ql) , (5) 

where fjp/p is the so-called Pomeron flux and /j/jp is the parton distribution function of the 
Pomeron. 

Note that the QCD fits to the diffractive data lead to aip(O) for the effective Pomeron 
trajectory, which is somewhat larger than the one for the effective soft Pomeron trajectory. 
This is likely due to a different interplay of soft and semihard physics in hard diffraction 
at the Qq scale and a different role of screening compared to soft interactions. Hence, it is 
likely that a violation of the factorization approximation will be be observed once the data 
are more accurate. 

The final expression for f^/^^^ takes the form 



fj'/fi^^ Ql, ^p) = -^^^'^Bjf]p/p{xip) fj/]p{(3 = x/xip, Ql) J d% 
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One immediately sees from Eq. (jH)) that the factorization approximation is not valid for 
nuclear diffractive parton distributions, even if it is valid for the nucleon case: At fixed xp, 
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the right hand side of Eq. © depends not only on P but also on the Bjorken x since the 
screening factor is given by the exponential factor containing alg, which is a function of x. 
In addition, the right hand side of Eq. ^ depends on the atomic mass number A since the 
effect of nuclear shadowing increases with increasing A. 

The aforementioned breakdown of the factorization approximation is a result of the in- 
crease of the nuclear shadowing effects both with the increase of incident energy and with 
the increase of the atomic number. This precludes the possibility of a scenario offered in 
Ref. PH], where coherent diffraction in DIS on nucleon and nuclear targets is provided by 
the same universal diffractive PDFs -"a universal Pomeron". 

It is also worth noting that the approximation which we use at Qq in order to take into 
account multiple rescatterings, corresponds essentially to treating diffraction as superposition 
of elastic scattering of different components of the photon wave function off the nucleus. This 
is a reasonable approximation for the configurations with masses comparable to Q^. As one 
approaches the /3 -C 1 limit (which corresponds to Mx ^ Q^, one approaches the limit 
analogous to the soft triple Pomeron limit, in which case diffraction off nuclei is strongly 
suppressed as compared to the elastic scattering, see e.g. [Sllin]- This effect should be even 
stronger in our case of DIS since increases with the decrease of (3. Hence, we somewhat 
overestimate diffraction at small (3 and at relatively small Qq scale, see Figs. |21 and El At 
larger Q^, diffraction at small (3 is dominated by the QCD evolution from /5 > 0.1 at Qo 
and, hence, the accuracy of our approximation improves. Hence, in the numerical studies, 
we neglect the effect of the small-/? suppression. 

3 Numerical results 

In our analysis of Eq. we used the 1994 HI fit PP for fj/ip, where the gluon distribution 
is decreased by the factor 0.75. This change seems to be required by the more recent analysis 
of the 1997 HI data on inclusive diffraction on hydrogen [2]. 

In order to have an idea about the magnitude of diffraction in DIS on hydrogen, the 
ratios Rj/n = ff/N^/fj/N for u-quarks and gluons and F^^'^^ / F2N are presented in Fig. [T]as 
functions of Bjorken x. Note that by definition 



These ratios give the probability of diffraction in the processes dominated by the coupling 
of a hard probe to a quark or gluon, respectively. For the hard process with a specific 
trigger, the probability of diffraction maybe close to Rq/N, such as the measurement of the 
diffractive structure function F.^, or to Rg/N, such as the 6-quark production. Alternatively, 
the probability of diffraction can have an intermediate value between Rq/N and Rg/N, such as 
in the s-quark production. In Fig. Q the solid curves correspond to Q = 2 GeV; the dashed 
curves correspond to Q = 10 GeV; the dot-dashed curves correspond to Q = 100 GeV. Since 
the dependence of Qn^"^^ I 9n is rather strong, we also show the ratio Qn^"^^ I Qn at Q = 5 GeV 
(dotted curve). One sees from Fig.^that in the quark channel, diffraction constitutes 15-20% 
of the total cross section, while in the gluon channel diffraction is significantly larger jlO] . 
This is likely to be related to a larger cross section of the interaction of the gluon color dipole 
(in the 8x8 representation) as compared to the triplet quark- ant iquark dipole. 




(7) 
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Figure 1: The ratios ff/^/fj/N for the M-quarks and gluons and NLO F^^'^^/F2n- The 
sohd curves correspond to Q = 2 GeV; the dashed curves correspond to Q = 10 GeV; the 
dot-dashed curves correspond to Q = 100 GeV. In addition, for the gluons the dotted curve 
correspond to Q = 5 GeV. 



The absolute upper limit of the gluon distribution, /qn = 1/2 is reached at a; = 
6 X 10~^ and Q = 2 GeV. Since we prefer to stay away from modeling the kinematics, 
where taming of the increase of the diffractive parton distributions becomes necessary, we 
will consider the limited range of Bjorken x, x > 6 x 10~^, in this paper. 

Next it is natural to analyze how the ratios presented in Fig. ^ change when hydrogen 
is replaced by a nuclear target. This can be done in two steps. First, acting in the spirit 
of the QCD factorization theorem for hard diffraction in DIS, Eq. (jH)) is used to define the 
input for DGLAP evolution at fixed xp. Subsequent QCD evolution enables us to determine 
f^^^\x,Q'^,X]p) as a function of (3 = x/xp and at all fixed xp. An example of these 
results in presented in Figs. |21 and 121 for the nuclear targets of ^''Ca and ^°^Pb. The u- 
quark and gluon nuclear diffractive parton distributions f^/^^\x, Q^, xjp) (arbitrary absolute 
normalization) are presented as functions of /3 = x/xjp at two fixed values of xjp = 10~^ 
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and xjp = 10^^. The solid curves correspond to Q = 2 GeV; the dashed curves correspond 
to Q = 10 GeV; the dot-dashed curves correspond to Q = 100 GeV; the dotted curves 
correspond to Q = 5 GeV. Different shapes and sizes of f^^^\x,Q^,X]p) at xjp = 10~^ 
and xjp = 10~^ clearly demonstrate violation of the factorization approximation for nuclear 
DPDFs. 
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Figure 2: The w-quark and gluon nuclear (^°Ca) diffractive parton distribution as a function 
of f3 at two fixed values of xp. The solid curves correspond to Q = 2 GeV; the dashed curves 
correspond to Q = 10 GeV; the dot-dashed curves correspond to Q = 100 GeV. In addition, 
for the gluons the dotted curve correspond to Q = 5 GeV. 



Another characteristic feature of nuclear DPDFs is that, like in the case of the free proton 
target, the gluon distribution is significantly larger than the quark distribution. However, 
we point out that the ratio of the quark to the gluon DPDF is significantly larger in the 
nuclear case because of a faster increase of the quark nDPDF with the atomic number A. 
Also, similarly to the free proton case, scaling violations of nDPDFs at large (3 are rather 
insignificant. This point is exemplified in Fig. |3] where we plot the w-quark and gluon nDPDFs 
as functions of at fixed large (3 = 0.5 and small xjp = 10~^. One readily sees from Fig. |3] 
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Figure 3: The w-quark and gluon nuclear (^''^Pb) difFractive parton distribution as a function 
of (3 at two fixed values of xjp. The solid curves correspond to Q = 2 GeV; the dashed curves 
correspond to Q = 10 GeV; the dot-dashed curves correspond to Q = 100 GeV. In addition, 
for the gluons the dotted curve correspond to Q = 5 GeV. 



that QCD evolution in In is weak. 

The difference in large-/3 scaling violations of the structure function ^2^^^'* in the nuclear 
and nucleon case is presented in Fig. El At xjp = 10~^ and two values of j3, (3 = 0.5 and 
/3 = 0.1, ^2^^^^ is plotted function of for ^°Ca (solid curves), ^osp]^ (dashed curves) 
and free nucleon (dot-dashed curves). The curves are normalized to coincide at the lowest 
Q"^ = 4 GeV^. One can readily observe from Fig. El that scaling violations are largest for the 
free nucleon and that scaling violations decrease as one increases A. 

Having obtained f^^^^\x,Q'^,X]p), they can be integrated over over xjp at fixed Bjorken 

X, just like in Eq. ((7|). The resulting f^/^^ / fj/A ratios for the w-quarks and gluons and NLO 

F2A^^ / F2A are presented in Fig. (HI for ^°Ca and in Fig. [71 for ^°^Pb. From Figs. Inland [71 one 

can see that the fraction of the diffractive events in DIS at small x for moderately heavy 
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Figure 4: The u-quark and gluon nuclear PDFs as functions of Q"^ at f3 = 0.5 and small 
X]p = 10~^. The solid curves correspond to ^''Ca; the dotted curves correspond to ^°^Pb. 



and heavy nuclei is of the order of 30% and weakly changes with Q^, which is in a good 
agreement with the early estimates of Ref. jH]. In the case of gluon-induced reactions, the 
probability decreases rather significantly with an increase of Q^. However, the probability 
still remains at the level of 15-20% at Q = 10 GeV and, hence, it would be feasible to study 
this in ultraperipheral collisions at the LHC using, for instance, production of heavy flavors 
similarly to the case of inclusive production considered in [20]. Another option is to use dijet 
production like it was done in the proton case in the ZEUS [21j and HI j2] experiments. 

Another conclusion that can be drawn from Figs. El and [7| is that the A-dependence of 
the probability of coherent diffraction is rather weak for A > 40. For these values of A, the 
interaction for the central impact parameters is close to being completely absorptive (black) 
with a small contribution from the opaque nuclear edge. Moreover, the A-dependence is 
weaker in the gluon case since the gluon interactions at the Qo scale are already close to the 
black limit, even for the nucleon. 

Mathematically this pattern is a result of a compensation of two effects - stronger small- 
X nuclear shadowing in the case of coherent diffraction compared to the inclusive case, is 
compensated by the nuclear form factor, as a consequence of nuclear coherence. 
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Figure 5: The diffractive structure function F^''^^ as a function of at = 0.5 and (3 = 0.1 
and small xjp = 10^^. The solid curves correspond to "^^Ca; the dotted curves correspond to 
^°^Pb; the dot-dashed curves correspond to the free nucleon. 



It is worth noting a qualitative difference between the A-dependence of the fraction of the 
diffractive events in the quark and gluon-induced processes at small x. In the gluon case, it 
is a very weak function of A because already in the proton case, the probability of diffraction 
is close to one half, the maximal value allowed by unitarity. At the same time in the quark 
case a steady growth with A is predicted since for the proton the probability of diffraction 
in this channel is rather small and, hence, the increase of the blackness of the interaction 
with A leads to a gradual increase of the diffraction probability to the values close enough 
to the black body limit. 

4 Conclusions and Discussion 

We study small-x coherent diffraction in DIS on nuclear targets using the theory of lead- 
ing twist nuclear shadowing and the QCD factorization theorem for hard diffraction. It is 
demonstrated that Bjorken x and A-dependent nuclear shadowing explicitly breaks down 
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Figure 6: The ratios f^j^ I fjlA for the M-quarks and gluons and NLO F^^^ IF2A for ^°Ca. 
The sohd curves correspond to Q = 2 GeV; the dashed curves correspond to Q = 10 GeV; 
the dot-dashed curves correspond to Q = 100 GeV. In addition, for the gluons the dotted 
curve correspond to Q = 5 GeV. 



Regge factorization in diffraction, which means that at fixed xjp, nuclear parton distribution 
functions depend not only on /3 = x/xjp but also on Bjorken x and A. 

We calculate nuclear DPDFs (Figs. El El and Ej) as functions of /5, Q^, xjp and the atomic 
number A. Like in the free nucleon case, the gluon nDPDF is much larger than the quark 
nDPDF. Using the calculated nDPDFs, the contribution of coherent diffraction to the total 
probability is estimated for the w-quark and gluon channels as well as for the NLO F2 
structure functions (Figs. El and (71). The key result is an observation of dramatically different 
patterns of the A- dependence. In the quark channel and in the F^ case, the probability of 
diffraction increases with A, reaching about 30% at x ~ 10"'* and Qo ~ 4 GeV^ and for 
A ~ 200. In the gluon channel, the probability of diffraction is already large for the proton 
and, hence, it changes (decreases) rather insignificantly when the proton target is replaced 
by the heavy nuclear target: the probability remains at the level of ~ 40% at x ~ 10~^ and 
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Figure 7: The ratios f^j^ I fjlA for the u-quarks and gluons and NLO F^^^ IF2A for ^^'^Pb. 
The sohd curves correspond to Q = 2 GeV; the dashed curves correspond to Q = 10 GeV; 
the dot-dashed curves correspond to Q = 100 GeV. In addition, for the gluons the dotted 
curve correspond to Q = 5 GeV. 



^ 4 GeV2. 

The Q^-dependence of the probabihty of diffraction is also different in the quark and the 
gluon channels. The very large gluon diffractive distribution makes QCD evolution of the 
ratios u^'^^ juA and NLO F^a^ IF2A rather weak. At the same time, the ratio g^'^^ j Qa falls 
off rapidly as increases (compare the solid and broken curves in Figs. El and [7j). 

From the experimental point of view, coherent diffraction in deep inelastic scattering on 
nuclei can be identified via a two-step procedure. First, similarly to the case of ep scattering, 
one selects events with a rapidity gap. Second, one needs to separate the coherent and 
incoherent diffraction. This can be readily done using the lack of neutrons in the zero angle 
neutron calorimeter since the break-up of the nucleus in incoherent diffraction results in 
production of several evaporation neutrons, see discussion in jTHI. In addition, the ratio of 
incoherent diffraction to coherent diffraction is expected to be ~ 0.1 — 0.15 ^H]- Hence overall 
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in the collider kinematics the task of selecting the diffractive channel without break-up of 
the target appears to be much easier in the nucleus case than in the proton case. The t- 
dependence of coherent diffraction originates primarely from the factor (F4(t))^ where -Fa(^) 
is the nuclear form factor. Hence average t are small and it is hardly possible to measure the 
t-dependence of the diffractive amplitude for the case of the large masses of the produced 
diffractive system. However since the t-dependence is mostly trivial, inability to measure 
the differential cross section would not lead to a significant loss of information about the 
dynamics of diffraction. Note also that that the break-up channel originates mostly due to 
the scattering off the edge of the nucleus, leading to the same pattern of diffraction as in the 
scattering off a free nucleon. Hence we predict a different /3, xjp, dependence of the hard 
diffraction in incoherent and coherent diffraction. 

Nuclear diffractive PDFs, discussed in this paper, exhibit novel effects, which are large 
enough to be measured in the ultraperipheral collisions at the LHC. We also would like to 
emphasize that the proximity of the probability of hard diffraction to the unitarity limit at 

~ 4 GeV^ shows that the color transparency phenomenon and related to Bjorken scaling 
decomposition over powers of l/Q^, which are typical for DIS, disappear in the vicinity of 
these Q^. Thus, nuclear shadowing does not preclude observation of a variety of phenomena 
characteristic for the unitarity limit for the gluon channel in the case of the nucleon and 
nuclear targets at x < 10~^, and for the quark channel for x < 10~^ for heavy nuclear 
targets. In particular, we expect that blackening of the interaction will reveal itself in heavy 
ion collisions at the LHC (and to less extent at RHIC) in the filtering out of nonperturbative 
QCD effects and producing a pQCD phase in the proton-nucleus collisions in the proton 
fragmentation region [22\ and in the heavy ion collisions in the ion fragmentation regions [23] . 

Numerical results presented in this paper are available from V. Guzey (vadim.guzey@tp2.ruhr- 
uni-bochum.de) upon request. 

This work was supported by Sofia Kovalevskaya Program of the Alexander von Humboldt 
Foundation (Germany) and Department of Energy (USA) and GIF. 
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